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Abstract 

We present the first study of the statistics of GHz photons in quantum circuits, using Hanbury- 
Brown and Twiss correlations. The superpoissonian and poissonian photon statistics of thermal 
and coherent sources respectively made of a resistor and a radiofrequency generator are measured 
down to the quantum regime at milliKelvin temperatures. As photon correlations are linked to 
the second and fourth moments of current fluctuations, this experiment, which is based on current 
cryogenic electronics, may become a standard for probing electron/photon statistics in quantum 
conductors. 

PACS numbers: 73.23.-b,73.50.Td,42.50.-p,42.50.Ar 
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The seminal experiment of Hanbury-Brown and Twiss (HBT)^, consisted in two detec- 
tors correlating the power fluctuations of a single electromagnetic source at the outputs of a 
beam splitter (Fig^a). Positive correlations observed for thermal photons within the coher- 
ence time have been interpreted as originating from photon bunching due to Bose-Einstein 
statistics [2]. Since correlations are very sensitive to statistics, HBT experiments have been 
widely used to investigate optical sources: laser light has poissonian fluctuations revealed by 
vanishing HBT correlations; negative correlations have also been reported for non-classical 
sources j3|. 

HBT correlations have also been observed with degenerate electrons in quantum con- 
ductors, where the analogue of the photon flux is the electrical current. Due to Fermi- 
statistics, electrons emitted by contacts can be noiseless. As a result, sub-poissonian noise 
and electron antibunching are observed , in quantitative agreement with pre- 

dictions 

mm. 

Comparing photon and electron statistics would remain a formal 
exercise unless one notices the dual representation of the current as fermionic excitations in 
a quantum conductor and bosonic electromagnetic modes in the external measuring circuit. 
This rises the intriguing question on how to relate the sub-poissonian statistics of electrons 
fluctuations in the conductor to the statistics of the emitted photons QHQ- A flrst 
theoretical answer has been given in Ref. jl5l] where the photon statistics is shown to re- 
markably deviate from the superpoissonian black-body radiation. As shown below there is 
also a direct relationjDetween photon statistics and the fourth moment of the electron Full 
Counting Statistics jlJ ]. 

But, how to measure reliably the photon statistics in quantum circuits? Quantum effects 
in conductors, best displayed at sub-Kelvin temperatures, require few photon number de- 
tection at GHz frequencies and make experimental detection of TEM photon fluctuations 
challenging. An elegant approach is provided by Quantum Dots and Superconducting meso- 
scopic photon detectors ll^ • However, a more versatile approach is to use cryogenic low 
noise ampliflers (LNA) followed by detectors. 

In this letter, we present the flrst HBT photon correlation measurements in the GHz range 
and in the few photon population number regime. The use of linear phase- insensitive am- 
pliflcation prior to photon detection provides a new situation having no optical counterpart. 
It is shown to preserve the nature of the HBT correlations characterizing the photon source 
under test. The sensitivity to various statistics is tested. First, an impedance matched resis- 
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tor is used which provides an ideal black-body thermal source. Although expected, we show 
for the first time that the Johnson-Nyquist noise is associated with superpoissonian noise 
and positive HBT correlations characteristic of Bose-Einstein thermal distribution. Second, 
a low noise radiofrequency generator is used as a monochromatic coherent photon source. It 
is found to generate coherent photons with poissonian statistics and vanishing HBT corre- 
lations (like a Laser). Regarding auto-correlations, which are sensitive to detection details, 
a quantum description of the amplifier noise is made which quantitatively agrees with our 
measurements. 

Following Nyquist the current noise power of a conductor of conductance G can be 
related to the power of TEM waves emitted in the external circuit. If Z is the characteristic 
circuit impedance, the average power in a small frequency range z/, v + dv is: 

^ Z 



[l + GZ] 



;Si{i')di' = Nhudu 



Here is the photon population, and Si{i>) = {^ly/du the spectral density of current 
fluctuation ^I{t) filtered in the frequency range v + dv The fluctuations of the 

current noise power entail fluctuations in the photon power at frequencies much lower than 
V. For classical currents, its variance, proportional to the low frequency bandwidth is: 



Gzy 

and that of the photon population is directly related to the fourth moment of the current 
fluctuation: 



Ze' V (A/)2 



^h{l + GZyj {e^uduy 
According to Glauber [2i|, thermal fluctuations of a classical current lead to superpoissonian 



fluctuations (AA^)^ = A^(l + N) or: 

((AP)2) = 2BN{1 + N){hufdu (1) 

with A^ given by the Bose-Einstein distribution. In the case of quantum conductors, the 
current is no longer a classical observable and the unsvmmetrized current noise operator 
Si{h',t) = dr i{t)i(t + T)exp{i2'7ih'T) is needed |l^|2l|. A full quantum treatment can 
be found in Ref. which shows that, for conductors with large numbers of electronic modes, 
the statistics of photons shows only small deviations from Bose-Einstein distribution even 



when shot noise dominates over thermal noise. However, a single electronic mode conductor 
with transmission close to 1/2 is found to emit non-classical photons for frequency in the 
range I/e. Non-classical photons are characterized by negative deviations (ANy — from 
the Poisson statistics, which is precisely the quantity that HBT correlations measure. This 
remarkable result and the connection with the fourth moment of the full electron current 
statistics provide a strong motivation to investigate HBT photon correlations in quantum 
circuits. 

The radio-frequency equivalent of the optical HBT experiment is shown in Fig^b. Source 
(a) is the rf-photon source under test. The emitted TEM photons propagate through a short 
50 Ohm characteristic impedance coaxial line and are fed to a cryogenic 3dB strip-line power 
splitter. The splitter scattering matrix, measured with a network analyzer, is identical to 
that of an optical separatrix (phases included). A built-in 50 Ohms resistor, source (b), 
plays the role of the vacuum channel of the optical case. Photon vacuum is achieved when 
its temperature Tq satisfies To ^ Tq where Tq = hu/kB- Outputs (1) and (2) of the power 
splitter are not immediately detected but a 1-2 GHz linear phase-insensitive amplification 
chain is inserted before square-law detection (see Fig^c). Each chain consists in a mi- 
crowave circulator followed by an ultra-low noise cryogenic amplifier and room temperature 
amplifiers. The circulators, at low temperature, ensure that amplifiers do not send back pho- 
tons towards (a) and (b). The detectors give an output voltage proportional to the photon 
intensity after amplification P°'^. Their finite 1 /is integration time allows to monitor the 
low frequency photon intensity fluctuations. A fast numerical spectrum analyser calculates 
the autocorrelations ((AP°2*)^) ^-^d the HBT cross-correlations (AP°"*AP2°"*) in the band 
40-200 kHz {B =160 kHz). 

In a first series of experiments, A and B, the sensitivity to Bose-Einstein statistics is 
tested. In a third experiment C, Poisson's statistics is tested using coherent photons. In 
A, source (a) is a 50 conductor whose temperature T is varied from 20 mK to several K, 
while source (b) realizes good photon vacuum (Tq = 17 — 20mK< Tq). A pair of 1.64- 
1.81 GHz filters select a narrow band frequency around v =1.72 GHz, with Tq = 86 mK 
and hdu <^ ksT. A quantum amplifier description (see below) predicts the mean powers 
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Pi = P"^^/Gi, referred to the input, and their fluctuations: 




t = 



1, 2 



(2) 



2 



(3) 



Gi 2 2 hi^ 



(AP.AP,) = 2BC-f - fnhuy 



(4) 



Here Na and Nb are the photon populations of sources (a) and (b) given by Bose-Einstein 
distribution, with A^;, as To << Tg. Autocorrelations differ from Eq.(^ in two ways: 
flrst, the factor 1 in the second parenthesis, revealing the independent particle behavior of 
thermal photons is replaced by 1/Gi {Gi = 80 dB) when refered to the amplifler input; sec- 
ondly, an extra photon population is added due to amplifler noise expressed in temperature 
units Tjv,i (Tjv ~ 15K in the experiment A and T/v ~ 6K in experiments B and C). The 
HBT cross-correlation are unaffected by ampliflcation. 

Experimental results are shown in Fig|21 The solid line is the Bose-Einstein theoretical flt 
of the mean photon power with Eq.Q taking T/v and G as free parameters. The experimen- 
tal data reproduce well the quantum crossover at Tq/2 ^43 mK. This is the coldest quantum 
cross-over ever reported for microwave photons. The parameter G is consistent with inde- 
pendent set-up calibration. The experimental scatter, 6T ^ 1 mK in temperature units, 
corresponds to the resolution 6T ^ 2T/v / VduAt expected for few seconds acquisition time 
At. This unprecedented sensitivity ammounts to flux rate variations of 1 photon per /is. The 
HBT cross-correlations (AP1AP2) are expected to vary like [Na) ~ (T/Tq) for T > Tq 
as Nb = 0. This is exactly what we observed in FiglHl This provides the flrst evidence for 
Bose-Einstein correlations of photon emitted by a resistor in the few photon number limit 
at sub-Kelvin temperature. The correlation resolution is 5 ((AP)^) ^ ((APi 2)^) / V BAt. 
This corresponds to detect population fluctuations of \fAN1KN2 ~ ksT^/ (^v\/ BAt^ ^ 
1.5 for our experimental At =1000 s. 

Experiment B provides evidence for vanishing HBT correlations (AP1AP2) ~ (^Na — iV^)^ 
when Na = Nb. The experiment is performed at higher temperature (T = 4-24 K, Tq = 4 
K). Source (b) is no longer in the ground state: Nb{i') ~ fc^To/Ziz/ ~ 40. Here, the central 
frequency is z/ = 1.5 GHz, the bandwidth du ^ 0.8 GHz and the amplifler noise temperatures 
TN,i = 6 K and 8 K. As shown in FigOl Pi,2 linearly depend on temperature in agreement 



5 



with Eq.(j2I). The T = 0-extrapolates define the LNA noise temperatures. In our analysis 
the hnear slopes give d =1.53 10^ (Az/i = 0.75 GHz), = 1.17 10« (Az/2 = 0.88 GHz) 
in fair agreement (within 20 %) with independent calibration. The square root of the auto- 
correlation shows linear variation with temperature exemplifying the super-poissonian noise 
of thermal photons given by EqI21 There are no free parameter left for the HBT correlations 
as gains are known. They are found positive and quantitatively agree with Eq. (jU (solid line 
in FigO)). They vanish for equal temperature sources. 

Experiment C tests the sensitivity to a different statistics. A microwave source of fre- 
quency 1^0 = 1.5 GHz and bandwidth 100 Hz generates monochromatic photons. Its output 
is attenuated at cryogenic temperature Tatt- The average power Pa in (a), is chosen compa- 
rable to that delivered by the thermal source in experiment B. As for a Laser, the source is 
expected to generate coherent photon states |2^. The photon statistics being Poissonian, 
the low frequency power fluctuations of source (a) in bandwidth B are ((APq)^) = 2Bhi'QPa. 
An important question is whether attenuation and arnpliflcation before detection change the 
statistics and if yes how? According to references 22,123], the commutation rules for bosonic 
input and output operators imply the addition of an extra bosonic operator describing the 
amplifler/attenuator noise. Applying this quantum constraint, we flnd (in units referred to 
the amplifier chain input): 

((APi,2)'> = 2Bhv^F^^J\2 (AP1AP2) = 
Pi,2 = 1 + 2kB r 

where we have used Tatt = Tq. Again, amplification has no effect on HBT cross-correlations 
and the absence of cross correlation characterizing Poisson's statistics remains. In the auto- 
correlations, a Fano factor F appears due to amplification noise. The results of experiment 
C are shown on FigEJ Cross- and auto-correlations are plotted versus detected power re- 
ferred to the input. Residual cross correlations, due to to small set-up imperfections, remain 
negligible at the scale of the autocorrelations. Auto-correlations are perfectly linear with 
power as expected. We measure a large Fano factor Pi^2 = 310, 400 in accordance with the 
order of magnitude of (Tatt + Tni^2) /hvo. 

In conclusion, the highly sensitive HBT photon correlations performed in the GHz range 
at subKelvin temperature using phase insensitive LNA and square-law detection easily dis- 
criminate between different statistics. In particular the cross-correlations are unaffected by 



amplification details. The method, simple, versatile and based on currently available elec- 
tronics, can be easily reproduced in other laboratories. It appears to be very suitable to 
study the photon population statistics of TEM modes emitted by quantum conductors or 
equivalently the fourth moment of current fluctuations. 
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Figure 1: Schematics of the HBT experiment (a), of the GHz beam spHtter (b), of the ampUficatioii 
and detection chains (c). A detailed description is given in the text. 
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Figure 2: Thermal source at u = 1.7 GHz. Main frame: quantum crossover at Tq/2 ^ 40 mK in 
the mean photon power Pi as function of temperature T, or effective occupation number kBT/hu 
for a beam sphtter temperature Tq 17 mK. The sohd hne is a Bose-Einstein fit with an amphfier 
noise temperature T^v ~ 15 K. Inset: positive HBT correlations show the bunching of thermal 
photons and the absence of residual correlations at T = K. Experimental scatter lies within the 
expected statistical uncertainty (error bar). 
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Figure 3: : Superpoissonian correlations of a thermal source in the classical regime: mean photon 
power Pi and noise (APjAPj) (bandwidth B = 160 kHz) as function source temperature T with 
a beam splitter at Tq « 4 K. The temperature scale of ((APj)^) and Pj is used as a calibration 
of the measuring line. Cross-correlations ^yJAPlA^j are deduced without extra parameter and 
found to vanish as (T — Tq) in agreement with Eq.Q. 
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Figure 4: Poissonian correlations of a coherent monochromatic source. Autocorrelations are linear 
in the input RF power with a large Fano factor F due to amplification and attenuation. Cross 
correlations vanish in accordance with theoretical predictions. 
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